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Abstract 



We present asymmetries between the production of D + and D~ 
mesons in Fermilab experiment E791 as a function of xf and pf. The 
data used here consist of 74,000 fully-reconstructed charmed mesons 
produced by a 500 GeV/c ir~ beam on C and Pt foils. The measure- 
ments are compared to results of models which predict differences be- 
tween the production of heavy-quark mesons that have a light quark 
in common with the beam (leading particles) and those that do not 
(non-leading particles). While the default models do not agree with 
our data, we can reach agreement with one of them, PYTHIA, by 
making a limited number of changes to parameters used. 
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Among the least understood aspects of quantum chromodynamics (QCD) is 
fragmentation, the non-perturbative dressing of bare quarks into hadrons. Certain 
asymmetries in the production of particles and antiparticles are especially sensitive 
to fragmentation effects, while largely free of experimental bias in the measurement 
process. In this paper we present high statistics measurements of charmed particle 
asymmetries from 7r-nucleon interactions at experiment E791 at Fermilab. 

Early studies of charm particle production |l|-|3| showed evidence of a large 
enhancement in the forward production of charmed particles that contain a quark 
or antiquark in common with the beam (leading particles) over those that do not 
(non- leading particles). Leading-order QCD calculations predict no asymmetry 
between quark and antiquark production. Several other models were developed 
to explain this effect although none fully accounted for the difference in 

production between leading and non-leading particles. As sample sizes and signal- 
to-background ratios increased, experiments ||[7| with signals on the order of a 
thousand charmed mesons continued to find a forward production asymmetry, one 
much higher than predicted even by next-to-leading-order (NLO) QCD calcula- 
tions which became available [f|. Now, we report results based on 74,000 fully 
reconstructed D ± decays, enough to examine the asymmetry both as a function of 
Feynman-x {xp = p z /pT ax ) an d transverse momentum squared (pf). 

For a 7T~ beam, NLO QCD calculations predict a small enhancement in the 
number of c quarks over c quarks in the very forward direction. However, as 
has been previously shown |7j] , this does not fully account for the observed leading 
particle enhancement. We compare our results to two current models which predict 
a much larger enhancement than NLO QCD. We also report changes to parameter 
values within one of these models (PYTHIA) |9|,[lIJ that provide agreement with 
our data. 

One possible explanation of the asymmetry is given by the beam-dragging 
model in which forward momentum is added to the produced heavy quark if it 
combines with a remnant light quark from the incoming beam particle, forming a 
leading particle. This causes leading charmed particles to have a harder Xp spec- 
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trum than non-leading particles. This feature is characteristic of the Lund string- 
fragmentation model which is the basis of the fragmentation simulation software 
JETSET used in PYTHIA. 

A second class of models [|lT|-|l^1 postulates charm may be produced as virtual 
cc pairs in the beam particle in addition to being produced perturbatively. An 
example of the former process is the intrinsic charm model proposed by Vogt and 
Brodsky fi~l|| . Here, the incoming 7r~ beam particle fluctuates into a \udcc ) Fock 
state. Since the virtual cc quarks have nearly the same velocity as the original u 
and d quarks in the pion, they are likely to coalesce with these remnant quarks, 
forming leading particles. Charmed quarks can also hadronize via fragmentation 
mechanisms which do not distinguish between leading and non-leading charm pro- 
duction. In reference jll] specifically, the fragmentation mechanism assumes that 
the momentum of the D meson is equal to the momentum of the c quark. 

To quantify the difference in the production of leading and non-leading parti- 
cles, an asymmetry parameter A is defined for each region of phase space: 

Mx»p!) = 2^01. (i) 

Ol + cr NL 



Here, <7 L (a NL ) is the production cross section for the leading (non-leading) particle 
being studied. Experiments WA82 |J and E769 |7[] found large values of A (~ 0.5) 
for high values of xp (— 0.6). 

For the ir~(du) beam in E791, the D°(cu) and D~(cd) are leading charmed 
mesons. However, a large fraction of the D°'s are produced by the D* + — > D°tt + 
decay process. The original D* + is actually a non-leading particle. Therefore, the 
observed D 0, s come from a mixture of leading and non-leading processes, making 
their study more complex. In order to minimize systematic errors, we restrict 
ourselves to studies of the asymmetries in the D ± spectra. 

The data were recorded from 500 GeV/c n~ interactions in five thin foils (one 
platinum, four diamond) separated by gaps of 1.34 to 1.39 cm. The E791 spec- 
trometer is an upgraded version of the apparatus used in Fermilab experiments 
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E516, E691, and E769 [14]. Momentum analysis is provided by two dipole mag- 
nets which bend particles in the horizontal plane and 35 planes of drift chambers. 
Silicon microstrip detectors (6 in the beam, 17 downstream of the target) provide 
precision track and vertex reconstruction. 

E791 recorded 2 x 10 10 events with a loose transverse energy trigger. After 
reconstruction, events with evidence of multiple vertices are kept for further anal- 
ysis. We search for D mesons in the decay mode D + (D~) — > K~ -k + ii + (K + n~ n~) . 
The selection criteria listed below are chosen to maximize Ng / \j ~Ns + Nb where 
Ns and Nb are the number of signal and background events. 

To help ensure that the reconstructed secondary vertex is a true decay vertex, 
we require that it be separated from the primary vertex by at least 9.0 07, where 
o\ is our resolution on the separation, and separated from the nearest target foil 
by at least 0.75 <j\. In addition, we impose two requirements to ensure that the 
candidate decay tracks do not originate at the primary vertex. First, each of the 
three decay tracks must have an impact parameter with respect to the primary 
vertex greater than 5.5 a p , where a p is the resolution on the impact parameter. 
Second, we form the ratio of each track's distance from the secondary vertex to 
that from the primary vertex and require the product of these three ratios to be 
less than 0.002. Another two requirements ensure that the reconstructed charged 
D is consistent with originating at the primary vertex. First, the impact parameter 
of the reconstructed D momentum vector with respect to the primary vertex must 
be less than 50 \mi ||15|| . Second, the component of the reconstructed D momen- 
tum perpendicular to the D line-of-flight (as determined from the primary and 
secondary vertex positions) must be less than 0.7 GeV/c. Finally, to help ensure 
that the three measured tracks come from the decay of a high mass particle, the 
sum of the squares of the tracks' momenta transverse to the direction of the recon- 
structed candidate momentum must be greater than 0.4 (GeV/c) 2 . The resulting 
Kttti mass spectra are shown in Fig. [I] for D + and D~ candidates. 

The Kirn spectrum, in each region of phase space, was fit to a Gaussian curve 
for the signal plus a linear distribution for the background. The number of D can- 
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didates is determined from the total area under the Gaussian curve. The width of 
the Gaussian was fixed for each region, with the evolution of the width as a function 
of xp determined by Monte Carlo simulations. The systematic error in the num- 
ber of mesons caused by the uncertainty in the width was estimated by increasing 
and decreasing the width by 10%, recalculating the yield each time. The resulting 
systematic uncertainties are small compared to the statistical uncertainties. 

Detailed studies of our charged particle reconstruction using both data and 
Monte Carlo simulations show a small difference in reconstruction efficiency be- 
tween positively and negatively charged particles due mainly to degradation over 
time of the wire tracking chambers in the vicinity of the incident beam. Because 
the negative beam particles are bent by the analysis magnets, the center of degra- 
dation was shifted slightly from the center of our detector. The effect on the 
measured values of A was determined using Monte Carlo simulations and is signif- 
icant only in the high xp region, increasing the value of A by ~ 0.15 in the highest 
bin. The systematic uncertainty caused by the efficiency correction was estimated 
by modeling the degradation of our chambers for four different time periods during 
our run, including periods near the beginning and the end of our data-taking. The 
rms spread of the asymmetries from the four different efficiency corrections was 
used to estimate the systematic uncertainties. These uncertainties are comparable 
to the uncertainties resulting from statistical errors. 

Figure |2| shows the value of the asymmetry parameter calculated for different 
bins of xf , compared to results from experiments E769 and WA82. The E791 
error bars shown on all plots include both systematic errors (from the uncertainty 
on the mass resolution and from the uncertainty on the reconstruction efficiency) 
and statistical errors, added in quadrature. Our measured asymmetries are higher 
than in those two experiments, especially at low xf, perhaps due to the different 
beam momenta (500 GeV/c for E791, 340 GeV/c for WA82 and 250 GeV/c for 
E769). 

In Fig. | we compare our values of A to predictions based on two models for D 
production previously discussed and to the prediction for charm quarks by NLO 
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QCD. The results from the Lund string fragmentation model used in PYTHIA 
are shown in two different curves. The dashed line shows the prediction from the 
default PYTHIA and is significantly higher than our data for —0.2 < Xp < 0.4. 
This is due to the fact that PYTHIA predicts a higher overall production ratio 
of D~ to D + than is seen in data. We also explored a range of parameters to 
determine if our results could be accommodated within the PYTHIA model. A 
tuned PYTHIA prediction []IB| that is consistent with our data is shown as the 
solid line in Fig. |3|. This tuning included increasing the mass of the c quark (m c ) 
from 1.35 GeV/c 2 to 1.7 GeV/c 2 and increasing the average primordial k\ of the 
partons from (0.44 GeV/c) 2 to (1.0 GeV/c) 2 . This decreases the likelihood that the 
remnant d quark can combine with the c quark with a small enough invariant mass 
to collapse to a D" meson. While this value of mean /c 2 may seem unphysically 
high, similar values of k\ are suggested by other observations from photo- and 
hadro-production charm experiments [17| . 



In Fig. [3] we also compare our results to a recent prediction involving intrinsic 
charm by Vogt and Brodsky, specifically calculated for a 500 GeV/ c n~ beam 
[|ITyi"8|l . Although the shape of the A vs xp curve is similar to our data, the 
prediction is too low for all xp. This model assumes equal numbers of D~ and D + 
mesons were produced, unlike the PYTHIA model, and changing this ratio may 
produce better agreement with data. 

Figure f| shows our data compared to the predictions for A from these same 
models versus pi for —0.2 < xp < 0.8. Again, the default PYTHIA model predicts 
values which are too high for most points, while the tuned PYTHIA is in good 
agreement with our data. The intrinsic charm model of Vogt and Brodsky predicts 
values close to zero, indicative of the assumption that D~ and D + mesons are 
produced in equal numbers. 

In the intrinsic charm model, coalescence occurs dominantly at low pi (~ m 2 ) 



lTf . Thus, in the context of coalescence alone, A should be large at high xp and 
low p 2 , increasing as p 2 goes to zero. In order to address this issue, we show the 
asymmetry as a function of pf for mesons with high xp values (0.4 < xp < 0.8) in 
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Fig. |5|. Here, the intrinsic charm model predicts values of A that are too low to 
match our data. However, we note that fragmentation mechanisms alternative to 
the one used in reference [11] can soften equally the xp spectrum of both the leading 
and non-leading Z)'s. This would emphasize the hardness of the Xp spectrum of 
leading Z)'s produced by coalescence and thereby increase the predicted asymmetry 
at high xp. Independent of this over all scale, the model predicts an increase in 
the value of A as pf goes to zero. Our data show no indication of this increase. 

In summary, we see a substantial leading-particle effect at high xp which cannot 
be accounted for using current perturbative QCD calculations. The beam-dragging 
model used in PYTHIA can account for the asymmetries between D~ and D + 
production if reasonable changes are made to the parameters used in the model. 
Neither the default PYTHIA nor the intrinsic charm model is consistent with our 
data. 
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FIG. 1. Ktttt mass distributions for D + and D~ candidates after all selection criteria 
are applied. 
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FIG. 2. The charged D production asymmetry A as a function of xf- The E791 
results are compared to previous results from WA82 [6] and E769 [7]. WA82 data are 
for D ± from a 340 GeV/c n~ beam, E769 data are for D ± and D*± from a 250 GeV/c 
7r ± beam, and E791 data are for D ± from a 500 GeV/c ir~ beam. E791 and E769 data 
are for mesons with pi less than 10 (GeV/c) 2 . 
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■ E791 Data 
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NLO QCD (Frixione) 




J I I I I I l_l I I I I I I I l_l I I I I I I I I I I I I I I I I I L_l 



J I I I I I I I L 



-0.2 -0.1 







0.1 



0.2 



0.3 0.4 0.5 0.6 0.7 0.8 

X F 

FIG. 3. The asymmetry A as a function of xf for E791 data (points with 
error bars), for three models, and for NLO QCD. The tuned PYTHIA prediction [16] is 
from the default PYTHIA Monte Carlo software with the c quark mass increased to 1.7 
GeV/c 2 and the average primordial k\ of partons increased to (1.0 GeV/c) 2 , as discussed 
in the text. The Vogt-Brodsky prediction is specifically calculated for the E791 beam 
momentum [18]. Both data and model predictions are for mesons with pi less than 
10 (GeV/c) 2 . The NLO QCD prediction is for charm quarks rather than D mesons [8]. 
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FIG. 4. The asymmetry A as a function of pf for our full range of xf 
(—0.2 < xp < 0.8). The data points correspond to the measured values. The curves 
represent theoretical predictions. The tuned PYTHIA prediction [16] is from the default 
PYTHIA Monte Carlo software with the c quark mass increased to 1.7 GeV/c 2 and the 
average primordial k\ of partons increased to (1.0 GeV/c) 2 as discussed in the text. The 
intrinsic charm prediction is specifically calculated for the E791 beam momentum [18]. 
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FIG. 5. The asymmetry A as a function of pf for high values of xf 
(0.4 < xp < 0.8). The data points correspond to the measured values. The curves 
represent theoretical predictions. The rise in A at low values of pf predicted by the 
Vogt-Brodsky model is a key feature which is not reproduced by our data. 
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